Capillary-driven flow induced by a stepped perturbation atop a viscous film 
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Thin viscous liquid films driven by capillarity are well described in the lubrication theory through the thin 
film equation. In this article, we present an analytical solution of this equation for a particular initial profile: 
a stepped perturbation. This initial condition allows a linearization of the problem making it amenable to 
Fourier analysis. The solution is obtained and characterized. As for a temperature step in the heat equation, 
self-similarity of the first kind of the full evolution is demonstrated and a long-term expression for the excess 
free energy is derived. In addition, hydrodynamical fields are described. The solution is then compared to 
experimental profiles from a model system: a polystyrene nanostep above the glass transition temperature 
which flows due to capillarity. The excellent agreement enables a precise measurement of the capillary velocity 
for this polymeric liquid, without involving any numerical simulation. More generally, as these results hold 
for any viscous system driven by capillarity, the present solution may provide a useful tool in hydrodynamics 
of thin viscous films. 



INTRODUCTION 

Micro- and nanofilms are of tremendous importance in 
a variety of scientific fieldJ^l, such as polymer physics, 
physiology, biophysics, micro-electronics, surface chem- 
istry, thermodynamics or hydrodynamics. For instance, 
they are involved in modern mechanical and optical engi- 
neering processes, through lubrication, paints and coat- 
ing. Gaining a complete understanding of these sys- 
tems is a key step towards the development of molecular 
electronics, biomimetics, superadhesion and self-cleaning 
surfaces. 

Thin films furthermore remain of fundamental inter- 
est in physics and mathematics, as they raise important 
questions that are still unsolved. The example of polymer 
systems, which we will preferentially refer to throughout 
the present article, is enlightening on this very point. As 
far as rheology of ultra-thin polymer films is concerned, 
when the height of the film is comparable to the char- 
acteristic size of the macromolecule, scaling arguments 
have been proposed for the effective viscosity 4 . Sev- 
eral of the recent experimental results that have been 
obtained motivate the desire for a fundamental under- 
standing of polymers in confinement. Examples are the 
enhancement of the effective mobility in the liquicP and 
glassy state^, and long-term relaxation of bulk viscositjd. 
The modification of p olym eric conformations^ and in- 
terchain entanglements^^ near surfaces, and in con- 
fined geometries^, have been studied in detail. Sur- 
face instabilities and pattern formation have been ex- 
plored as welf^HHI More over, the role of film prepara- 
tion has been invest igat e d 1 5 * 1 6 1 but remains complex, as 
pointed out in th e particular case of the glass transition 
temperatur e 1 17 * 1 ^ . 
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FIG. 1. A schematic showing the capillary leveling of an 
initially stepped perturbation with amplitude So atop a thin 
viscous film of height ho. The vertical profile h(x,t) of the 
free surface depends only on the horizontal coordinate x and 
time t. 



The systems discussed above are often well described 
by the lubric ation theory through capillary-driven thin 
film equations^!. However, due to their high orders 
and non-linearities, these particular equations have not 
yet been solved analytically. Nevertheless, the reader 
will find details on the mathematical advances in recent 
reviews^20|. In addition, we point out the fact that thin 
film equation s hav e been solved numerically in various 
configurations^^. 

In the present communication, we address the analyt- 
ical problem of a stepped perturbation atop a flat film 
(shown schematically in Fig. [I]) which is infinite in the 
two horizontal dimensions. The interest in this system 
is twofold. First, it is physically relevant since it is in- 
spired by experiment! 23 * 24 ! . Secondly, it is mathematically 
interesting since the smallness of the step in comparison 
to the underlying film height allows one to linearize the 
thin film equation and solve it through Fourier analysis, 
as for the heat equation 2 ^. 

This study is divided into three sections. In the first 
one, we recall the main ingredients of the physical model, 
including the general non-linear thin film equation for 
two-dimensional capillary-driven flows. In the second 
part, we explicitly consider the stepped perturbation, for 
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which we obtain and characterize the linear solution. Fi- 
nally, in the third part, we compare the theoretical results 
to experiments performed on a polystyrene (PS) stepped 
film above the glass transition temperature. 



the case of inhomogeneous surface tension is considered 
elsewhere^' and that we neglect confinement effects^^ 
since the film height is still large in comparison with the 
size of the macromolecule, typically ~ 10 nm for the con- 
sidered molecular weight. 



I. PHYSICAL FRAMEWORK 

In this first section, we describe the physical model 
that is considered. After justifying the main assump- 
tions and giving the boundary conditions for the flow, 
we derive the general non-linear thin film equation in 
two dimensions. 



A. Assumptions 



We consider the example of PS nanofilms above 
their glass transition temperature T g ~ 100 °C, since 
this experimental system will be compared to the the- 
ory in section |III| We thus have the following typical 



B. Vertical boundary conditions 

We consider the case of no shear at the liquid-air in- 
terface: 



d z v\ z=h = 



(4) 



where z is the vertical coordinate. In addition, we assume 
no slip at the substrate: 



V\\\z=0 = , 



(5) 



with uii the component of the velocity parallel to the 
liquid-substrate interface. 



parametergSSHSST vertical height ho ~ 1 um, dynam- 
ical viscosity rj ~ 1 MPa.s, molecular weight M w ~ 

15 kg.mol" 1 , surface tension 7 ~ 30 mN.nT 1 , density C Thin film et l uation 

p ~ 1 g.cm~ 3 and shear modulus G ~ 1 MPa. Extension 
to any thin viscous fluid is straightforward by using the 
corresponding relevant orders of magnitude. 

Let us first estimate the typical spreading capillary ve- 
locity v c , Reynolds number Re, capillary length l c and 
Maxwell viscoelastic time tm- 
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From these orders of magnitude, and since we typi- 
cally o bserve slow evolution of the profiles over several 
minute^ 23 ! 24 !, we can make the following assumptions: we 
have an incompressible flow of a highly viscous Newto- 
nian fluid, where gravity^, disjoining pressure®! and in- 
ertia are negligible. Therefore, this flow is well described 
by the Stokes equation: 

VP = ?)Av , (2) 

combined with the incompressibility condition: 



In addition to the assumptions and the vertical bound- 
ary conditions presented in the previous parts, we assume 
a spatial invariance along one horizontal direction y. We 
then have a pure two-dimensional problem. Therefore, 
the height of the free surface is given by h(x,t), where x 
is the relevant horizontal direction and t the time. The 
local pressure field is a priori given by P(x, z, t). Accord- 
ing to the lubrication approximation, we can neglect the 
vertical velocity with respect to the horizontal one and 
write: v = v{x,z,t) e x , where e x is the horizontal basis 
vector. Then, we project and integrate Eq. pj), using 
Eq. p}, Eq. Q and Eq. (j5j, and we find: 



d z P = 



(6) 



the local pressure field P(x, t) is thus invariant in the 
vertical direction, and: 



v(x,z,t) = ^-(z 2 



2hz) d x P , 



(7) 



which corresponds to the usual parabolic Poiseuille flow. 
Volume conservation requires that: 



d t h + d x I dz v = 
Jo 



(8) 







(3) 



where P and v are the local pressure and velocity fields 
in the liquid. In addition, we assume that the profile 
slopes remain small in comparison to 1, which is an in- 
gredient of the lubrication approximation. Finally, we 
take 7 and rj as homogeneous and constant. Note that 



Finally, because the pressure does not depend on z (see 
Eq. (|6|) it can be evaluated at the free surface through 
the Young-Laplace equation. Since the lubrication ap- 
proximation implies small curvatures, the pressure satis- 
fies: 



P-Po 



-7<9 2 /i , 



(9) 
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where Pq is the atmospheric pressure. Thus, combining 
Eq. 0, Eq. ^ and Eq. ^, we get: 



d t h+^d x (h 3 d*h) =0 , 



(10) 



which is knowrP ^ * 15 ! 21 * 31 ! as the capillary-driven thin film 
equation. 

II. STEPPED PERTURBATION 

In the present section, we consider an infinitesimal per- 
turbation 6(x, t) <C ho of the free surface on an infinitely 
large and flat thin film with height h (see Fig. [T]). We 
thus linearize the thin film equation, before nondimen- 
sionalizing the problem and giving its general formal so- 
lution through Fourier analysis. Then, we consider a sim- 
ple initial profile: a stepped film, for which we express 
and characterize the solution. In particular, we demon- 
strate self-similarity of the first kind of the full evolution 
and we study the long-term viscous dissipation of the ex- 
cess free energy, before describing the hydrodynamical 
fields of the problem. 



A. Linearized Thin Film Equation 

In the following, we restrict ourselves to the linearized 
version of the thin film equation. Setting: 



h(x, t) = ho + S(x, t) , 



Eq. ( 10 ) becomes: 



d t 6 + ^d*6 = 0, 
3rj 



(11) 



(12) 



to first order in the perturbation. Apart from the higher 
order in the spatial derivative, this equation is analogous 
to the heat equation^. 



C. Formal general solution 

Injecting a propagating mode e l ( KX ~ nT ) ; f angular 
frequency Q and angular wavenumber K, into Eq. ( 14 ) 
leads to the dispersion equation: 



(15) 



which implies in particular that spatial oscillations with 
small wavelengths decay in amplitude faster than those 
with large wavelengths. Since Eq. ( 14 1 is linear, we 



can write the general solution as a superposition in the 
Fourier basis. Using Eq. (151, it follows that: 



A{X,T) = I -^L A(K) e- K ' 1 e 1 



—KTAKX 



(16) 



where A(K) is the spatial Fourier transform of the initial 
profile: 



A{K) = I -= A(X,0) e" 1 



KX 



2tt 



(17) 



D. Solution for a stepped initial condition 

We consider now the particular ideal case where the 
initial profile is a step proportional to the Heaviside func- 
tion 0: 

A(Jf,0) = A e(X) , (18) 

through a dimensionless amplitude A = Sq/Hq. 

Under the initial condition of Eq. (JTsJ) , Eq. (16 1 be- 
comes: 



A(X,T) = ^[1 + V(X,T)] 



(19) 



where we shifted the vertical origin and normalized the 
amplitude through the centered profile: 



B. Non-dimensionalizing 



<P(X,T) 



KX 



(20) 



For generality, we introduce the typical height ho and 
time t = 3r)h /j of the problem, as well as the associated 
dimensionless variables: 



and where the dashed integral represents Cauchy's prin- 
cipal value. 





5 


A = 






ho 




X 


X = 






ho 




t 


T = 






to 



We non-dimensionalize Eq. ( 10 1 and obtain: 



<9 T A + <9 X 4 A = 



(13a) 
(13b) 
(13c) 

(14) 



Limits and symmetry 



First, Eq. (20) naturally reduces to the sign function 
at T = 0, as required from the initial condition through 
Eq. (18). Secondly, the energy cost associated with the 



surface perturbation is viscously dissipated and the final 
equilibrium state is thus flat at finite X: 



lim A(X,T) 



Ao 
2 



(21) 
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Thirdly, at finite T, using Cauchy's residue theorem and 
Jordan's lemma, and provided we can perform an inver- 
sion of limits, we obtain: 



lim A(X,T) = A Q 

X— 5- + 00 

lim A(X,T) = 0. 



x 



(22a) 
(22b) 



This important result tells us that the horizontal bound- 
ary limits remain those of the initial profile all along the 
evolution, which stresses the fact that the boundary lim- 
its of the initial profile are crucial in defining the sym- 
metry of the shape. Moreover, the leveling does not con- 
taminate the infinite horizontal limits at finite time, as 
expected. 

Finally, there is a fixed point at X = 0, along the 
temporal evolution: 



A(0,T) 



Ao 
2 



(23) 



which is a center of symmetry because ip is an odd func- 
tion with respect to X. 



Self-similarity 

At finite T, let us change variables through: 
X = UT 1/4 

Q 



K = 



JU/4 ■ 



Thus, Eq. (20) becomes: 



1 p(X,T)= X (U) 



+ 00 



dQ 



3 -<2 4 
inQ 



(24a) 
(24b) 

(25a) 
(25b) 



Therefore, with the initial profile considered, the full evo- 
lution is self-similar of the first kincP^: it depends only 
on the variable U and the profile satisfies the following 
equatiorPS 



X 



U 



x 



(26) 



whose i ntegra l has been studied in the mathematical 
literature^!]. Once again, this equation is analogous to 
the one obtained with self-similar solutions in XT^ 1 ! 2 for 
the heat equation^, apart from the order of the deriva- 
tive on the left-hand side. Note that the self-similarity 



introduced in Eq. ( 24 ) and Eq. ( 25a ) is highly dependent 
on Eq. ( 22 I, and thus on the boundary limits of the initial 
profile. Different boundary limits would lead to different 
forms of self-similarity^. In fact, we see through Eq. ( 16 ) 



that the condition for such a self-similarity of the solu- 
tion is that the Fourier transform of the initial profile is 
positive homogeneous of degree —1: 



A(aK) = oT x A{K) 



(27) 



for a > 0. In other words, according to Eq. (17) 
A(aX,0) = A(A,0) . 



(28) 



Therefore, the initial profile has to be a constant func- 
tion with a potential discontinuity at X = 0, i.e. the 
self-similarity considered can only be achieved if there is 
a step initially, which justifies a posteriori the interest of 
this configuration. In addition, note that these results are 
analogous to the ones obtained for the heat equation^! 
with a temperature step, except for the precise time expo- 
nent which is related to the order of the spatial derivative 
in the governing equation. In this picture, we see that 
the final equilibrium state (see Eq. ( pT| )) and fixed point 
(see Eq. (23)) values are unified into a single point of the 



self-similar profile: 

X(0) - . (29) 

Moreover, the profile is symmetric with respect to this 
particular central point: 



x(-u) = - x (u) 



(30) 



Finally, we see through Eq. ( 25b ) that \ is proportional 
to the inverse Fourier transform of the function Q t— > 
Q-i e -Q _ Therefore, the solution is: 



(31) 

1Z7T 



where we introduced two auxiliary functions for clarity: 



F a (17) = lF 3 



F (U) = X F 3 
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4J 
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4J 





(32a) 
(32b) 



with the de finitio n of the (1, 3)-generalized hypergeomet- 
ric function^ 37 ! 38 ! 



1 F 3 ({a} : {b,c,d},w) = J2 



(a) k 



fe>0 



(b)k(c)k(d)k k\ ' 



(33) 



and the Pochhammer symbol of the rising facto- 
rial. Note that the self-similar dimensionless solution in 



Eq. (31 ) is unique: it depends only on the variable U and 
thus not on the experimental parameters ho, So, 7, 77 and 
t. 

Fig. [2] shows A/A as a function of X for various 
times. As we can see, the profile flattens through time as 
the excess surface energy is viscously dissipated. Apart 
from the oscillatory behavior linked with the fourth spa- 
tial derivative of Eq. (14), this evolution is qualitatively 



close to the solution of the heat equatiorPSl for which the 
same analytical treatment would lead to the usual error 
function. At fixed time T, the spatial oscillatory behav- 
ior (see inset of Fig. [2]) is qualitatively analogous to the 
temporal response of a damped harmonic oscillator after 
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FIG. 2. Solution of Eq. 



14) for a stepped initial profile 

to 



18)) at various dimensionless times 



(see Eq 

Eq. (Iffy , Eq. (24a} , Eq. (25a} and Eq. {51]). The profiles are 
symmetric with respect to the central fixed point, according 
to Eq. \3(ty . The self-similar profile is naturally identical to 
the real profile at T = 1, since x(U) = tp(U, 1) according to 
Eq. \24a\l and Eq. , 25a). The inset shows a zoom of the first 



oscillations for X > 0, in the T — 1 profile. 



having switched on an external constant load, although 
the spectrum is different according to Eq. ( 20 1 . Finally, 



plotting the same data with respect to the self-similar 
variable U, instead of variable X for a given parameter 
T, would allow to collapse all the data onto the A(Z7, 1) 
profile since x(U) = tp(U,l), by construction. 



G. Extrema 

Let us now consider the amplitude of the oscillatory 
behavior. A local extremum at U* of the self-similar 
profile is defined by x'{U*) — 0. So, we write: 



dQ e 



-Q 4 e iQ(7* 



. 



(34) 



It is straightforward to see that x'(~U*) = as well, and 
that x"{—U*) — —x"{U*). Therefore, each local max- 
imum (minimum) has a mirror minimum (maximum), 
whose horizontal coordinate is symmetric with respect 
to U = 0. This is trivial when recalling the symmetry of 
the profile through Eq. (30). Interestingly, Eq. (34) im- 
plies that U* is independent of Ao: the amplitude of the 
perturbation only shifts the vertical origin and stretches 
vertically the self-similar function x- This scaling in- 
variance of the solution is not expected to hold for the 
general non-linear thin film equation, due to the cubic 
dependency on the profile height (see Eq. (10)). 



We now focus on two mirror extrema indexed by i and 
whose coordinates are: [U*,x(U*)} and [— U* , x(— U*)], 



with U* > 0. By introducing the ratio: 
x (U*)-x(+oo) 



n 



x(-Ut)-xi-oo) 



(35) 



and by invoking the symmetry of the profile through 
Eq. (30), we naturally get: 



K = 1 



(36) 



which is indeed observed in Fig. [2] An important feature 
is that this ratio does not depend on the amplitude of 
the perturbation A . This scaling invariance of the so- 
lution does not hold for the general non-linear thin film 
equation^, due to the cubic dependency on the profile 



height (see Eq. (10)) 



H. Excess free energy 

In order to characterize the relaxation to equilibrium, 
we study the long-term evolution of the capillary free 
energy (per unit length along y) in excess of the equi- 
librium flat profile. At large time and thus small slopes, 
this quantity satisfies: 



A J" : 



\ I dx {djf 



(37) 



Rewriting Eq. (37) in dimensionless variables through 
Eq. (pi) gives: 



A.F 



jho 
2 



dX (d x Ay 



(38) 



Finally, invoking the variable U from Eq. (24a) leads to 
the desired long-term expression: 



AT 7 _ t(T\h ' 



(39) 



where we introduced the initial excess capillary energy 
AJ^ = "/So, and the characteristic time: 



as well as a constant numerical factor: 

31/4 



dU x (U) « 0.16 



(40) 



(41) 



The final relation co mes f rom a calculation of the integral 
J x' 2 ~ 1: using Eq. (31). Although the geometrical de- 
pendencies of this result are only limits of the non-linear 
case, the dependencies on material properties do hold in 
the non-linear case^. This important long-term result 
teaches us that the surface energy is viscously dissipated 
with a 1/4 power-law in the inverse time, and that the 
initial geometry plays the crucial role of a driving ampli- 
tude. 
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I. Hydrodynamical fields 



Using the solution obtained in Eq. (31 1, we now de- 



rive the expressions of the pressure and velocity fields in 
real variables, as well as the corresponding volumetric 
flow rate. As for the previous treatment of energy, using 
0, Eq. Jjj}, Eq. JIT), Eq. p|, Eq. p|, Eq. p4a" 



Eq. 

and Eq. ( |25a ), we find 



P(x,t)=P -x"(U) 



1/2 



and: 



3 3/4 

^*,M)« — x "'(t/) 



7*i 



1/4 



2z 
/it, 



"0 



(42) 



(43) 



to lowest order in the perturbation, 
velocity at the central fixed point x - 



In particular, the 
satisfies: 



v (z,t) 



-0.1 



7* 



1/4 



2z 
ho 



"0 



(44) 



where we calculated the numerical factor through 
Eq. (31 1. The minus sign is expected due to the chosen 



geometry (see Fig.JT]). Interestingly, we see that the am- 
plitude of the parabolic velocity profile, i.e. the surface 
velocity, decreases with time, viscosity and film height, 
and that it increases with surface tension and amplitude 
of the perturbation. Then, by integrating Eq. (43) over 



z, we obtain the volumetric flow rate (per unit length in 
the y-direction) : 



i x >"(u) W 1 

2 3V4 I v t 3 



(45) 



to lowest order in the perturbation. Thus, at the central 
fixed point x = 0, we have: 



Qo(t) 



-0.07 



rjt 3 



1/4 



The dependencies are the same as for the velocity field, 
except for the film height: when ho tends to infinity the 
flow diverges, as expected, whereas the velocity tends to 
zero. 




FIG. 3. Comparison between theory (see Eq. (3D) and ex- 
periments (see Eq. The experimental profiles corre- 
spond to AFM lines taken at room temperature after quench- 
ing from 120° C. PS stepped films were prepared as detailed 
elsewher P 3 ^ 24 \ The experimental data represent steps with 
different heights {ho, 80} in nm, as indicated. The profiles 
are recorded after t — 10 min of annealing; exceptions are 
for: {1063, 100}, where there are three profiles corresponding 
to t = 6, 8 and 10 min; and for {924,57} where there are 
two profiles corresponding to t = 10 and 20 min. There is 
only one free horizontal stretching parameter in this compar- 
ison, which leads to a measurement of the c apill ary velocity: 
7/77 = 1.5 ±0.1 um.min - , according to Eq. (51). 



A. Self-similarity 

Knowing the profile 8(x,t) and the film height ho for 
a given experimental sequence, it is straightforward to 
obtain the experimental equivalents (in real variables) of 
Eq. (pal): 



A C (M) 



( 46 ) and of Eq. (p|: 



5(x,t) 



(47) 



(48) 



We checked the self-similarity by plotting the experimen- 
tal profiles as a function of the real self-similar variable^ 



III. COMPARISON WITH EXPERIMENTS 



tV4 



(49) 



In this last section, we compare the analyti cal so lution 
obtained in section [ll| to a model experiment! 23 * 24 ! : a PS 
stepped perturbation above the glass transition temper- 
ature (see Fig. [TV In particular, we fit the profiles to the 
solution of Eq. pT| . As we finally see, this procedure en- 
ables an accurate measurement of the capillary velocity 
of the material. 



When doing so, the profiles at different times collapse 
onto a single profile: 



(50) 



thus demonstrating the self-similarity of the first kind of 
the experimental evolution. This result appears in Fig. [3] 
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B. Fitting the profiles 

When comparing the profiles x e ( u ) horn the exper- 
iment (see Eq. (50)) and x(U) from the theory (see 
Eq. (31)), we immediately see that a horizontal stretch 
factor remains due to the nondimcnsionalizing procedure. 
In fact, according to Eq. (13), Eq. (24a) and Eq. (49), we 
have: 



U 



3i] 



1/4 



(51) 



We thus fit all the experimental profiles x e ( u ) to the the- 
oretical one x(U) through this single free parameter. The 
result is plotted in Fig. [3] for various ho, So and t. 

The agreement is excellent and holds for height ratios 
up to So /ho ~ 10%. In particular, we see that the ex- 
perimental profiles plotted this way do not depend on 
ho, S and t, as suggested by the theory. The fact that 
the scaled profiles are independent of the stepped film 
geometry confirms the uniqueness and self-similarity of 
the dimcnsionlcss solution. Furthermore, since we know 
ho experimentally, this single parameter fitting procedure 
offers a precise measurement of the capillary velocity 7/77 
of the material at the considered temperature. We find: 
7/7/ = 1.5±0.1 um.min - 1 at 1 20°C, which compares well 
with the tabulated values^^U through the WLF modeP^l. 

As a final remark, we note that the initial stepped pro- 
file seems in contradiction with the lubrication approxi- 
mation used in the theory. Thus, the short-term evolu- 
tion of the experimental profile may not be described by 
the theoretical solution in Eq. (31 1 . However, as observed 



experimentally, the step levels due to Laplace pressure 
gradients and reaches small slopes on a timescale that 
is short in comparison to the one of all the experimen- 
tal observations. Furthermore, the experimental profiles 
connect rapidly to the theoretical solution. Therefore, 
the possible temporal offset is negligible in comparison 
with the typical experimental times and does not reduce 
the precision on the capillary velocity measurement. 



CONCLUSION 

We reported on an analytical solution of the linear thin 
film equation for a stepped initial condition. The solution 
was obtained by performing Fourier analysis and involved 
generalized hypergeometric functions. We characterized 
the solution and demonstrated the self-similarity of the 
first kind of the full evolution. In addition, we found 
the origin of this symmetry in the initial profile itself. 
Then, using self-similarity, we derived the long-term vis- 
cous dissipation law in such a linear system. This scal- 
ing captures the relevant physical ingredients and pro- 
vides a limit for the corresponding expression in the non- 
linear theory. We also derived the scaling expressions for 
the pressure and velocity fields and for the volumetric 
flow rate. As far as the linearized equation is concerned, 



apart from the well known case of the heat equation, we 
suspect that hypergeometric solutions may be obtained 
in an identical way for higher even orders of the spa- 
tial derivative in the linear equation. However, the odd 
orders being free from dissipation are expected to lead 
to fundamentally different mathematical solutions. Fi- 
nally, we compared the solution to experimental profiles 
obtained on polystyrene stepped films for several times, 
and height ratios up to ~ 10%. The agreement is excel- 
lent, thus demonstrating the interest of such a solution 
in the physics of capillary-driven thin viscous films. In 
particular, the fitting technique offers a precise viscome- 
ter without invoking any numerical simulation. In the 
near future, these results may be extended to other initial 
conditions and to the first non-linear term of the equa- 
tion through perturbation theory. The goals would be 
to study the convergence to intermediate self-similarity 
and to capture the non-linear extension of the present 
results in order to approach the full non-linear solution, 
which is still an open analytical problem that governs a 
considerable number of exciting physical applications. 
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